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I. Introduction 

The precise tests for the electroweak scale physics have established plenty of 
knowledge about the elementary particles [H [2J. The standard model (SM) has 
been evidenced to be indeed a very successful theory at the electroweak energy 
scale. Nevertheless, there are a lot of the issues in the flavor physics and universe 
observations for which the SM is not able to account [3]. During the past more 
decades a series of experiment results of B physics and neutrino physics have told 
us a great deal of the detailed information of flavor physics [I]. What are paid 
more attention are some facts in the following. Firstly, the fermion mass spectrum 
emerges a very large hierarchy. The quark and charged lepton masses range from 
one MeV to one hundred GeV or so [TJ, while the neutrino masses are only at the 
Sub-eV level [5]. Secondly, the quark flavor mixing pattern is distinctly different 
from the lepton one. The former are only three small mixing angles [6], whereas the 
latter has bi- large mixing angles and a bit small but non-zero # 13 [7]. In addition, the 
CP violation in the quark sector has been verified to be non-zero but very small [S], 
however, it is in suspense whether the CP violation in the lepton sector vanishes 
or not [9]. Thirdly, that the light neutrino nature is Majorana or Dirac particle 
has to be further identified by the experiments such as 0vj3j3 [10J. On the other 
hand, the astrophysics observations and researches lead to some impressive puzzles 
in the universe, in particular, the genesis of the matter-antimatter asymmetry and 
the original nature of cold dark matter [TTJ. Finally, a very important and unsolved 
problem is whether flavor physics are truly connected with the baryon asymmetry 
and/or cold dark matter in the universe or not? All the problems have a great 
significance for both particle physics and the early universe [12], so they always 
attract great attention in the experiment and theory fields. 

The researches for the above-mentioned problems have motivated many new 
theories beyond the SM. The various theoretical approaches and models have been 
proposed to solve the intractable issues [13J. For instance, the Froggatt-Nielsen 
mechanism with the U(l) family symmetry can account for mass hierarchy [14], 
the discrete family group A4 can lead to the tri-bimaximal mixing structure of the 
lepton mixing matrix [T5] , the non-Abelian continuous group SU(3) is introduced 
to explain the neutrino mixing [12], the model with the SO (3) family group can 
accommodate the experimental data of the quarks and leptons by the fewer param- 
eters [17] . some grand unification models based on the SO(10) symmetry group can 
also give some reasonable interpretations for fermion masses and flavor mixings [18] . 
In addition, some suggestions of the baryon asymmetry and cold dark matter are 
very constructive [19] . Although these theories are successful in explaining some 
specific aspects of the above problems, it seems very difficult for them to solve the 
whole problems all together. Indeed it is especially hard for a model construction 
to keep the principle of the simplicity, economy and the less number of parameters, 
otherwise the theory will be excessive complexity and incredible. On all accounts, a 
better theory beyond the SM has to be confronted with the integration of particle 
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physics and the early universe, in other words, it should simultaneously account for 
flavor physics, the baryon asymmetry and cold dark matter. Of course, it is still a 
large challenge for theoretical particle physicists to find a desirable theory to uncover 
these mysteries of the universe [20J . 

In this work, I construct a practical model to integrate three parts of the fermion 
flavor, cold dark matter and leptogenesis. The model has the local gauge symme- 
tries SU(3) C ® SU(2) L ® U(l) x <8> U(1) B -l- The subgroups U(l) x ® U(1) B -l only 
appear above the TeV scale, and then are spontaneously broken to the supercharge 
subgroup U(1)y below that scale. Besides the SM particles, some new non-SM par- 
ticles are introduced into the model. They are the fourth generation quarks and 
leptons, two vector gauge fields related to U(l) x ® U(1)b-l, and three scalar fields, 
which are respectively one neutral singlet, one charged singlet and one symmetric 
triplet under SU(2) L . The U(l)x <S> U(1)b-l breaking is accomplished by the neu- 
tral singlet scalar developing a non-vanishing vacuum expectation value (VEV) at 
the TeV scale. This breaking leads that one of the two new gauge fields obtains a 
TeV scale masse by the Higgs mechanism, in the meantime, the fourth generation 
neutrino is given rise to a TeV scale mass. The fourth generation neutrino has some 
unique natures in the model, which ensure that it is a cold dark matter particle. 
All of the SM particle masses are generated after the electroweak breaking. The 
triplet scalar field takes part in the symmetry breakings but it develops only a tiny 
VEV by virtue of its very heavy mass. This is a real source of the tiny masses of 
the light neutrinos. The leptogenesis is implemented by the charged singlet scalar 
decaying into a SM charged lepton and a cold dark matter neutrino. The decay 
process simultaneously satisfies the three conditions of the B — L violation, CP 
violation and being out-of-equilibrium. This mechanism can naturally generate the 
B — L asymmetry, subsequently it is converted into the baryon asymmetry through 
the electroweak sphaleron process [ST]. In the model, the flavor physics is intimately 
associated with the cold dark matter and the matter-antimatter asymmetry in the 
universe. The model can completely accommodate and fit all the current experimen- 
tal data of the fermion masses and mixings, the cold dark matter and the baryon 
asymmetry, furthermore, it also predicts some interesting results. Finally, the model 
is feasible and promising to be tested in future experiments. I give some methods 
of searching the non-SM particles of the model in the experiments at the LHC [22J. 

The remainder of this paper is organized as follows. In Section II I outline the 
model. In Sec. Ill, the model symmetry breakings, the particle masses and mixings 
are introduced. In Sec. IV, the leptogenesis and cold dark matter are discussed. In 
Sec. V, I give the detailed numerical results. Sec. VI is devoted to conclusions. 

II. Model 

The local gauge symmetries of the model are the direct product groups of 
SU(3) C ® SU(2) L ® U(l) x ® U{1) B -l- The subgroup symmetry of U(l) x ® U(1) B -l 
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only appears above the TeV scale, and it will be broken to the supercharge subgroup 
U(1)y below that scale. The model particle contents and their gauge quantum num- 
bers are in detail listed in the following, 

G«(8, 1,0,0), Wj(l,3,0,0), ^(1,1,1,0), 1,0,1), 

Qi(3,2,0, -) , [u iR ,u 4R ,u 4L \(3,l,l,-) , [d iR ,d 4 R,d AL ](3,l,-l,-), 

^(1,2,0,-1), [e iR ,e AR ,e AL \(l, 1, -1, -1) , i/ 4L (l, 1, 1, -1) , 
#(1,2,1,0), 0(1,1,-2,2), S~(l, 1,-2,0), A(l,3,0,2). (1) 

All kinds of the notations are self-explanatory as usual. X M and are two new 
vector gauge fields related to the gauge subgroups U(l) x and U{1) B -l, respectively. 
The last two numbers in the brackets are exactly the U(l)x and U(1)b-l charges. 
The fermion subscript i indicates the first three generation fermions, in addition, the 
model newly includes the fourth generation quarks and leptons. It should also be 
noted that the fourth generation left-handed fermions are all singlets under SU(2)l. 
The right-handed neutrinos are absent in (1). The reason for this is that the neutral 
fermions un, are considered as Majorana particles, so the corresponding v R are 
not independent fields but rather determined by the relation v R = CUZ 7 ', in which 
C is a charge conjugation matrix. Besides the SM doublet Higgs field, two complex 
singlet scalar <fi and S~ , and a symmetric triplet scalar A are introduced in the 
model. After the gauge symmetry breakings, and S~ have respectively a zero 
electric charge and a negative electric charge. These scalar field representations are 
such as 




V2 ' V H 

In addition, H = IT2H* will also be referred. In brief, all the non-SM particles in 
(1) compose the new particle spectrum, and they play key roles in the new physics 
beyond the SM. In general, the new particle masses are the TeV scale or above it, 
so they should appear in the TeV scale circumstances, for example, in the early 
universe. 

Under the model gauge symmetries, the invariant Lagrangian of the model is 
composed of the three parts in the following. Firstly, the gauge kinetic energy terms 
are 

^Gauge = ^puregauge + if-fD^f + (£>" H)\D „if ) 

+ (D^(D^) + (D"S-y(DpS-) + Tr[(^A)t(^A)] , (3) 

where / stands for all kinds of the fermions in (1). The covariant derivative is 
defined as 

D, = d, + i (gsG^Y + 92W; T - + gxX^f- + g B _ L YjL^>j , (4) 
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where A a and r % are Gell-Mann and Pauli matrices, Qx and B — L are respectively 
the charge operators of U(l)x and U(1)b-l, and g 3 , g 2 , gx, 9b-l are four gauge 
coupling constants. After the gauge symmetry breakings, some gauge fields will 
generate their masses and mixings by the Higgs mechanism. 
Secondly, the model scalar potential is given by 

Vscalar = H + X H (H^H) 2 + fiffi (f> + X^^f 

+ M 2 S S + S- + \ S (S + S-) 2 + M A Tr[A f A] + A A (Tr[A f A]) 2 

+ A 1 (i/ t i/)(0 t 0) + (\ 2 (rfH) + A 3 (0 f 0)) rr[A+A] 

+ (X 4 (H^H) + A 5 (0 t 0) + A 6 Tr[AtA]) (S + S~) 

+ X 7 (H T AH^ + h.c.) , (5) 

where all kinds of the coupling parameters are self-explanatory. In general, the size 
of the coupling coefficients are [Ai, A 2 , . . . , A 7] ~ 0.1. For [Mj, M A ] > 0, M s and Ma 
are original masses of the particles S and A, respectively. Their values are M$ ~ 1 
TeV and Ma ~ 10 5 TeV in the model. For [^ 2 H ,^ < 0, both the singlet <p and 
the doublet H will develop non-zero VEVs. is responsible for the B — L breaking 
U(l)x ® U(1)b-l U(l)y, in succession, H completes the electroweak breaking 
SU{2)l ® U{1)y —> U(l) em . The triplet A also takes part in the two breakings 
because it is involved in the couplings of the last term in (5). However, it is about 
to be seen that A has only a very tiny VEV due to the heavy Ma, therefore A only 
plays a secondary role in these breakings but it plays a key role in the generation 
of the tiny neutrino masses. Lastly, the singlet S~ does not participate in any 
breakings because of its VEV vanishing. In comparison with the SM Higgs sector 
[23] . the phenomena of the model scalar sector are variety and more interesting. 
Thirdly, the model Yukawa couplings are written as 



^Yukawa — [yiJ^ ,U 4L M U4 J I y upu\ _j I I ^ 

+ ( yt'pr ir ) ( z ) 

+ { T mt) ( ) ( - ) 



■T 



J V 

+ yldULS-O u] u R + t/ 4 d mz5 + O dt t/ R + ytWES + O e h R + h.c. . (6) 

These Yukawa couplings have apparently uniform and regular frameworks. M/ 4 are 
exactly the masses of the fourth generation quarks and charged lepton, which are 
directly allowed by the model symmetries, however, Mf A are about few TeV based on 
the model theoretical consistency. The coupling coefficients, y{, y(, y[, y{, are chosen 
some real numbers since their complex phases can be absorbed by the redefined 
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fermion fields. Since every coupling coefficients scale the order of magnitude of itself 
term, they are arranged such a hierarchy as \y{\ < \y{\ < \y{\ < \y$\ ^ 1 in the 
model. The coupling matrices, J, Y* , F*, 0\ characterize the flavor mixings among 
the four generation fermions. At present we are lack of an underlying understand for 
the fermion flavor origin, however, it is believed that some flavor family symmetry 
is embedded in an underlying theory at a certain high-energy scale, but it is broken 
at the low-energy scale. The coupling matrices should imply some information of 
the flavor symmetry and its breaking. Therefore I suggest that the flavor structures 
of the coupling matrices have such a style as 



pf= u 




yf =u '' 



f=u, 








a .f 


( 


a f 


1 




-a f 


1 


/ 


\ 






c u 






02/ 




\ 


1 J 






(7) 



The flavor structures are both simple and reasonable, in particular, there are only few 
flavor parameters. The size of the flavor parameters are normally [|a/|, |ci/|, |c2/|] 
~ 0.1. The majority of their complex phases can be removed by the redefined 
fermion fields. The remaining complex phases will become the sources of the CP 
violations in the quark and lepton sectors. In (6), the y{l couplings have evidently a 
full flavor symmetry among the first three generation fermions but they are relatively 
smaller. The y[Y^ couplings only keep such a discrete symmetry S2 as fi ^ — ^3 
between the second and third generation fermions. The couplings between the first 
three and fourth generation fermions, y^F* , break the flavor symmetry but they 
are relatively bigger. Lastly, the couplings involving the charged scalar S ± , y{0^ , 
are the smallest ones. After the gauge symmetries are broken spontaneously, all 
kinds of the fermion masses, m U) d. e , rn u , M Vi , will be generated by the correspond- 
ing couplings and the VEVs of H, A, <ft, respectively. Finally, I in particular point 
out that the matrices y{l + y^Y^ can be diagonalized by the unitary matrix U as 
follows 



U^(y(l + y{Yf)U 



f 
Hi 



Un 




y/2a f y( 









Vi + V2a f y( 








y{ + 2y{ 



(8) 



The mixing angles of Uq are 612 = 623 
from the tri-bimaximal mixing matrix 1 2 



45°, 9 



13 — 

For \y( 



0° 

< 



It evidently distinguishes 
< \y(\, the first and 
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second eigenvalues are approximately the same size, and the third one is the biggest. 
This property of (8) plays a key role in the neutrino mass and mixing in the model. 

In summary, the above features of the particle contents and Lagrangian are very 
important not only for the particle masses and mixings, but also guarantee the cold 
dark matter and leptogenesis in the model. In the following sections of the paper, it 
is about to be seen that z/ 4i has unique natures and plays a special role in the model. 
It is actually a cold dark matter particle. The leptogenesis is really implemented 
by the decay S~ — > + VI. In a word, the above contents form the theoretical 
framework of the model. 



III. Symmetry Breakings and Particle Masses and Mixings 

The gauge symmetry breakings of the model go through two stages. The first 
step of the breakings is U(l)x <S> U(1)b-l — >• U(1)y, namely the B — L breaking, 
in succession, the second step is SU{2)l (g) U{1)y — > U(l) em , i.e. the electroweak 
breaking. The former is achieved by the real part component of developing a non- 
vanishing VEV at the TeV scale, while the latter is accomplished by the neutral 
component of H developing a non-vanishing VEV at the electroweak scale. In 
addition, the neutral component of A also develops a very tiny but non-zero VEV 
owing of the last term couplings in (5). The scalar vacuum structures and the 
conditions of the vacuum stabilization are easy derived from the scalar potential 
(5). The details are as follows 

<p^f±p^, H^ \ „o+„„ |, il V2, |, S-^S- 







(4*) 








W V2 V ^h-XI ' { ^ V2 V 4 ^-A? ' 

(A) = — = ~ X 7 V ^ V H „ -MV^H (s - ) = q (Q) 

V2 v / 2(2Af| + \ 2 v 2 H + A 3 vJ) ~ 2V2M 2 Aeff ' ^ 7 ' 

The stable conditions include [fi 2 H , fi 2 ^ Ai] < 0, [A#, A^, A5, Aa, A 2 , A 3 , A 4 , A 5 , A 6 , AX^Xh- 
X 2 ] > 0, and Mj 3> [v^v'jj]. In addition, |Ai| should be sufficient small so that 
is one order of magnitude larger than vh, in this way, this ensures that the B — L 
breaking precedes the electroweak breaking. M^e// in (9) is a effective mass of the 
A particle when the breakings are completed (see (10)). Provided that Ma ~ 10 5 
TeV, V(/, ~ 2.5 TeV, vh ~ 250 GeV, and A7 ~ 0.1, thus this naturally leads to 
t>A ~ 0.1 eV, consequently, gives the tiny neutrino masses. Thus it can be seen 
that the tiny nature of the neutrino masses essentially originates in the very heavy 
Ma in the model. In this sense, this is a new form of the seesaw mechanism [25J. 
Finally, the S~ field has a vanishing VEV, so it does not actually participate in the 
breakings. In short, all the conditions are not difficult to be satisfied so long as the 
parameters are chosen as some suitable values. 
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After the model gauge symmetry breakings are over, the following massive scalar 
bosons, H°, (ft , S ± , A [A , A ± , A ±=t ], appear in the scalar sector. Their masses and 
mixings are such as 



tan29h 



Xiv^vh 



X H v 2 H 



2 ' 



M HO,<t>o = (/V^l + X H v 2 H ) =F |A^| - \hV 2 h \ a/1 + tan 2 26 h , 
M 2 Seff = M 2 S + -{\ A v 2 H + X 5 v 2 + \ 6 v 2 A ), 

M 2 Aeff = Ml+ l -{\ 2 v 2 H + X 3 vl + 2\ A vl), (10) 

where 6 h is the mixing angle between H° and 0°. Provided that A^ ~ Xh > |Ai|, then 
tan26h < 0.1, thus the two neutral boson masses are approximately M H o « v^A^ v# 
and M^o w y^X^v^. In a similar way, there is also a very weak mixing between 
H° and A , or between 0° and A . However, these very weak mixings in the scalar 
sector can all be ignored throughout. M Se ff and M^e// are respectively the effective 
masses of the particles S ± and A. Provided that M s ^ v^, then M Se ff ~ v^. In 
view of Ma 3> t^, obviously, there is M^e// ~ M^. At present, M H o has been 
measured by the LHC [26J, it's value is 125 GeV. However, the model predicts that 
M^o and Ms are about 1 TeV, it are quite feasible to find the two bosons at the 
LHC, but the A particles are too heavy to be detected. 

In the gauge sector, the gauge symmetry breakings give rise to masses and mix- 
ings for some of the vector gauge bosons through the Higgs mechanism. The breaking 
procedure of U{\)x ® U{\)b-l — ► U(l)y is such as 

9x9b-l , a gx 
g x = = , tanO g = , 

9 2 x + 9%-l 9B—L 
Qy = Qx + (B - L), Q N = -tanOgQx + cot6 g (B - L) , 
Bp = cosQgXp + sinOgYp , Z'^ = -sindgX^ + cos9 g Yp , 

m Bm = o, ^ = i|g^)|^ = ^, (n) 

where gi,B^,Q Y are respectively the gauge coupling constant, gauge field, super- 
charge operator of U(1)y, Z' is an obtained mass neutral gauge field, and Qn is a 
new charge operator related to Z' . There are two gauge parameters g± and tan8 g in 
(11), however, g\ is not a free parameter but determined by the electroweak relation 
<7i = yjATta/cosOw, only the mixing angle tanO g is a free parameter. In addition, 
the last equation in (11) implies Mz> ^ 2giv ( f > , so Mz> should be few TeV or so. It 
should also be pointed out that the mixing between Z'^ and Z^, which is the SM 
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weak neutral gauge boson, is very small. Their mixing angle is given by 



„. sirv i 6 a cos6 a v'l, . . 

tan2ff ^ 2^4 ' • < 12 > 

Because of vjj/v^ ~ 1CT 2 and sin6 w ~ 0.5, this mixing can indeed be ignored. 

Below the electroweak scale, the Yukawa sector becomes clear and simple since 
the fourth generation heavy quarks and charged lepton have decoupled. After they 
are integrated out from (6), the effective Yukawa coupling matrices of the three 
generation quarks and charged lepton are given by 

/ \ 

Y e f = uAe = y(l + y f 2 Y f + (y() 2 F f ®F f \ ® = J b f b* f b f . (13) 

V b) I J 

Each term physical meaning is very explicit. According to the standard procedures, 
in (6) the symmetry breakings give rise to all kinds of the fermion masses as follows 



M f=u , d ,e = ~^feff , M v = -^(Vil + VW\ M UA = -^j= y U4 . (14) 



Obviously, there is the mass hierarchical relation M v ~ 0.01 eV Mf =u ^ e < M U4 ~ 
ITeV. In addition, the hierarchical coefficients, \y[\ \y[\ \y[\, will lead to the 
hierarchical masses of the three generation quarks and charged lepton. On the other 
hand, there is not such a term as F v <g> F v ^ in M u , so M v is distinguished from 
Mf =U; d, e ■ This is a primary source that the lepton flavor mixing is greatly different 
from the quark one. In short, the interesting features of the fermion mass matrices 
dominate the fermion masses and flavor mixings. 

Finally, the fermion mass eigenvalues and flavor mixing matrices are solved by 
diagonalizing the above mass matrices. The quark and charged lepton mass matri- 
ces Mf =Ut d t e are hermitian, while the light neutrino mass matrix M u is symmetry. 
Therefore, the mass matrix diagonalizations are accomplished as such 

U ] U M U U U = diag (m u , m c , m t ) , U\M d U d = diag (m d , m s , m b ) , 

UlM e U e = diag (m e , m M , m T ) , U\M u U* n = diag (m ni ,m n2 ,m n3 ) . (15) 

In the light of the characteristic structures of Mf =U:d:£ , the mass eigenvalues of 
the quarks and charged lepton are certainly some hierarchy, and the flavor mixing 
matrices Uf =u>dfi are all closed to an unit matrix. In contrast, an exact solution of 
the M v diagonalization can be given by use of (8) as 

U n = U , m 2 n2 -m 2 ni = 2V2yW 2 a„vl, m 2 n3 - m 2 n2 » 2(y" 2 v A ) 2 . (16) 

Obviously, U n is completely different from U f= u ,d,e , moreover, the two mass-squared 
differences can explain the neutrino data very well. The above results are convenient 
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for the following numerical analysis. The flavor mixing matrix in the quark sector 
and one in the lepton sector are respectively defined by [27J ES] 

Ul U d = Uckm , Ul U n = U PMNS ■ diag (e^\ e^, l) . (17) 

The two unitary matrices Uckm and U pmns are parameterized by the standard form 
in particle data group [TJ. fli, 02 are two Majorana phases in the lepton mixing. All 
kinds of the mixing angles and CP- violating phases can be calculated numerically. 
Finally, all of the results can be compared with the current and future experimental 
data. 

IV. Cold Dark Matter and Leptogenesis 

The model can naturally and elegantly account for the cold dark matter and lep- 
togenesis after the model symmetry breakings are completed. The fourth generation 
Majorana neutrino z/ 4 ^ own some unique properties. It has been seen from the model 
lagrangian that v±l has only the three types of the couplings, Val^u^lZ' ^^4l</>°, 
and ~eRV± L S~ . Provided that the mass order M V4 < M s < [M U4 , M d4 , M e4 ], the only 
decay channel of v^ L is z/ 4£ — > eiUjdk via an off-shell boson S + because there are 
some weak mixings between the fourth generation quarks and the first three gener- 
ation ones. If the coupling coefficients are \y{\ ^ 10~ 9 in the last line of (6), then 
the v^l decay width is estimated as ^ 1CT 44 GeV. In other words, the v^l lifetime 
is actually two orders of magnitude longer than the now age of universe, therefore 
it becomes a relatively stable particle in the universe. On the other hand, a pair of 
z/ 4 £ can annihilate into other particle pair. The annihilate processes are mediated 
by either the gauge boson Z' or the scalar boson 0°. Because the v± L mass is de- 
rived from the B — L breaking, M V4 should be around one TeV. Consequently, the 
Majorana neutrino z/ 4 ^ is genuinely a weak interactive massive particle (WIMP), of 
course, it also belongs to one of the fewer species of particles which can survive from 
the early universe to the now epoch. Therefore is a good candidate of the cold 
dark matter [29] . 

The annihilate channels of z/ 4 ^ have two ways. The principal annihilate process is 
that a pair of v^l annihilate into all kinds of the SM fermion pairs by the intermediate 
gauge boson Z'^ as shown in the figure (1). The other annihilate process is that 
two z/ 4 £ neutrinos annihilate into two Higgs bosons H° by the TeV scale boson 0° 
mediating. Because Ai is smaller, the cross section of the latter case is normally far 
smaller than one of the former except for some Breit-Wigner resonance points, so I 
can ignore the latter and only consider the former. The annihilate cross section of 
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Figure. 1. The diagram of the cold dark matter annihilating into all kinds of 
the SM fermion pairs by the gauge boson Z'^ at the TeV scale . 



the figure (1) is calculated as follows 



a(u iL +v«,-+f+f)- 2564{s _ Mh ^ + {Tz , Mz , n 2. \j j^mf R 

Rf = 3 (1 5 + ^ 2_ } 



4 



..tf iJVl,, r, IliflVl... , 

+ 2Q N (f L )Q N (f R )-L(l ^ + (Qjv(/ l ) - g^v^)) 2 -^(1 

s s Mt 



,2 \f . / 4m 2 



fL,fR,V4L V Z M 
9 

R' f = {Q 2 N (h) + QIUr)) (i - Tjf ) + QQN(h)QN(f R ) 



where s = 4M i 2 4 / (1 — 1> 2 ) is the squared center-of-mass energy, v is the velocity of v^l 
in the center-of-mass frame. The sum for fi, Jr count all kinds of the SM fermions 
who are permitted by kinematics. In fact, all of the relatively lighter m/ in (18) can 
been approximated to zero except for the relatively heavier mt,M VA . On the basis 
of WIMP, the relic abundance of u^i in the current universe is determined by the 
annihilation cross section as such 

n ^ 2 58xr '' &V ", (19) 
(crv r ) 

where v r is the relative velocity of the two annihilate particles. In addition, the 
heat average (19) can be calculated by (crv r } rj a + b(v 2 } = a + bjp-, in which 
Tf w M U4 /20 is the freeze temperature of v^ L . A rough estimate is as follows. 
Because of g\ ~ 0.1, Yl Q%(f) ~ 10) ~ Mz' ~ 10 3 GeV, a weak cross section is 
naturally obtained as a ~ 10~ 9 GeV -2 , eventually, it leads to Qh 2 ~ 0.1, which is 
closed to the observation value. 
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Figure. 2. The tree and loop diagrams of the decay of S — > e i 
B — L asymmetry is generated . 



+ z/ 4 , by which the 



The baryon asymmetry through the leptogenesis can be implemented in the 
model. On the basis of the relevant couplings in (5) and (6), the main decay channel 
of the charged scalar boson S~ is S~ — > T^e^ (not including a heavy e~l ), in addition, 
S~ — > Widj can be ignored because its decay branching ratio is smaller. The figure 
(2) draws the tree and loop diagrams of S~ — > T>ie~[ , of course, its CP conjugate 
process S + — > has the corresponding diagrams. However, this decay process 
has simultaneously three items of the notable characteristics. Firstly, the decay 
evidently violates one unit of the B — L quantum number, namely A(B — L) = —1. 
Secondly, because there is a CP-violating source in the leptonic Yukawa sector, a 
CP asymmetry of the decay is surely generated by the interference between the tree 
diagram and the loop one. The CP asymmetry is defined and calculated as 



e 



T(S- -> e- + Vt) - T(S+ -» e+ + z/ 4 ) 
r[s- -> e; + vi) + T(S+ e+ + i/ 4 ) 

-v H X 4 £ £ M e Jm{OtT:pf){Ml - Ml - M^fj 

i=lj=l 

4V27rE|Of| 2 (Ml-il^ 4 ) 2 



a; 

~~ ± lv± v A 1V± S ~ lv± vi lv± s 



/j Xi - 1 ' ' MS. 1 M? - Ml. ' + Ml - Ml • { ' 



where provided with M e4 > Ms + M Ui . Thirdly, the decay is really an out-of- 
equilibrium process. Provided that the coupling coefficient \y%\ ^ 1CT 9 as before, 
then the decay rate is far smaller than the Hubble expansion rate of the universe, 
namely 



i-#EI?/^ 



16tt V M s 



H(T = M s ) hi 

M 



« 1 , (21) 
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where m 2 e jMg has been approximated to zero, M p i = 1.22 x 10 19 GeV, and g* is 
the effective number of relativistic degrees of freedom. In brief, the decay process 
S~ — > +V4 can indeed satisfy Sakharov's three conditions [30] . Therefore, this 
mechanism can naturally generate an asymmetry of B — L at the TeV scale. It can 
be seen from (20) that the CP asymmetry e linearly depends on the three quantities 
A 4 , Tij, fj, in particular, the major contribution comes from T i4 , namely the process 
by the ej inner mediating. In addition, e has no relation with y\ though the decay 
rate depends on y\. In the model, there are A 4 ~ 0.1, T i4 ~ 0.01, fj ~ 0.1, thus 
M s ~ 1 TeV will give e ~ 10~ 8 . 

The above B — L asymmetry occurs at the TeV scale. At that time the heavy 
S T and have been the non-relativistic particles, but the produced leptons are 
truly the relativistic states, moreover, their energy are normally more than 100 
GeV. Therefore the sphaleron electroweak transition can smoothly put into effect 
[31]. Consequently, the B — L asymmetry will eventually be converted into an 
asymmetry of the baryon number through the sphaleron process. According to the 
standard discussions, the baryon asymmetry is determined by 



where 7.04 is a ratio of the entropy density to the photon number density, c sp = 28/79 
is a coefficient of the sphaleron conversion, Yb-l stands for the B — L asymmetry, 
which is related to e by the expression in the parentheses, k is a dilution factor, 
it can actually be approximated to k w 1 on account of the very weak decay rate. 
At the TeV scale, only the SM particles are the relativistic statuses, whereas the 
non-SM particles are the non-relativistic statuses, so g* exactly equal to the effective 
number of degrees of freedom of the SM particles, namely g* = 106.75. In short, the 
baryon asymmetry can be calculated by the relations of (20) and (22). The model 
can achieve the leptogenesis at the TeV scale. 

V. Numerical Results 

In the section I present the model numerical results. The model involves a num- 
ber of the new parameters besides the SM ones. In the light of the III and IV section 
discussions, the parameters involved in the numerical calculations are collected to- 
gether in the following. The gauge sector has the gauge coupling gi and the mixing 
angle tan9 g . The scalar sector includes the three VEVs, v^, Vh, i>a> the two scalar 
boson masses Mho, Ms, and the scalar coupling A 4 . The Yukawa sector has M ei , 
all kinds of the coupling coefficients, and the flavor parameters, see (6) and (7). 
Among which, the three parameters, g±, vh, Mho, in fact belong to the SM parame- 
ters. Their values have been determined in the electroweak scale physics and by the 
recent LHC experiments [20], namely g\ = 0.356, vh = 246 GeV, M H o = 125 GeV. 
In addition, the three parameters, v^, va, M e4 , also belong to the fundamental pa- 
rameters in the model. Based on an overall consideration, the above six parameters 




(22) 
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are fixed throughout to the following values, 

0i = 0.356, v H = 246 GeV, M H o = 125 GeV, 

vj, = 2.5 TeV, v A = 0.1 eV, M e4 = 2.5 TeV. (23) 

The rest of the model parameters have to be determined by fitting all kinds of the 
experimental data. 

In the light of (6), the model Yukawa sector consists of the relatively independent 
quark sector and lepton one, so I will respectively discuss them. In the quark sector, 
the ten parameters in the following are involved in fitting the quark masses and 
mixing. Their input values are such as 

jtf = -1.1 x 10" 5 , V2 = -7.5 x 10~ 3 , 2/3 = 1, a u = -0.056, b u = 0.018 , 

y f = -l.l x 10~ 5 , ^ = 5.55xl0~ 4 , ^ = 0.153, a d = 0.168, b d = -0.039 i. 

(24) 

In the above the choices of the imaginary parts of a u , a d , b u , b d have a certain degree 
of freedoms. Here I only choose b d as a pure imaginary number since this style is 
simple and the best in the fits, of course, it is exactly the source of the CP violation 
in the quark sector. In short, this set of the values in (24) are reasonable, and 
also consistent with the previous discussions. According to the relevant relations in 
(13)-(17), the quark masses and mixing are solved as follows 

m u = 0.00234 , m c = 1.274 , m t = 173 , 
m d = 0.00481, m s = 0.0952, m 6 = 4.18, 

sf 2 = 0.2252, s| 3 = 0.0411, sf 3 = 0.00354 , 5 q = 0.377 tt = 67.8°, 

Jg P = 2.95 x 10" 5 , (25) 

where mass is in GeV unit, s a p = sin6 a p. In addition, the Jarlskog invariant J^ P , 
which measures the magnitude of the CP violation, is also figured out by using the 
quark mixing angles and complex phase. It is very clear that the numerical results 
in (25) accurately fit all the current experimental data of the quark masses, mixing, 
and CP violation [TJ. However, there are several points worthy pointed out in the 
successful fits. Firstly, the hierarchical parameters, y(, y(, y{, respectively dominate 
the 3rd, 2nd, 1st generation fermion masses, the parameters a d and a u impact on 
the mixing angles sin6i 2 and sin6i 3 , the parameters b d and b u are in charge of sin9 2 3 
and 5. Secondly, both the hierarchy of y(, y[, y( and the flavor structures of (7) 
jointly lead that the quark transform matrices U u and U d in (15) are both closed 
to an unit matrix, as a result, the quark mixing matrix XJqkm has eventually three 
small mixing angles. Thirdly, since the quark masses and mixing angles have been 
measured to a certain precision, the variable scope of the parameter space are very 
narrow. Fourthly, there is no any fine-tuning in the fits. Fifthly, there are the 
parameter relations yf = y", a d = —3a u , a d — a u = X c in (24), which A c = 0.224 
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is Cabibbo mixing angle. This should not be by coincidence, maybe there is an 
underlying reason for them. 

In the lepton sector, although the case of the charged lepton is similar to one of 
the quarks, the neutrino case makes great differences due to its particularity. It is in 
fact seen from (8) and (16) that U n has been fixed to Uo, and the three parameters, 
Hi, U2-, a u , determine the three mass eigenvalues of the neutrinos. Therefore, only 
y\, 1)2 are used to fit the two mass-squared differences of the neutrinos and a u can 
be fixed freely. The relevant parameters are chosen as follows 

^ = 3.51xl0~ 5 , ^ = 5.7xl0" 4 , yl = 0.0974, a e = 0.24, 6 e = —0.1 i , 
^ = 0.088, y£ = 0.304, a,, = 0.1, (26) 

where the pure imaginary b e is chosen as the source of the CP violation in the lepton 
sector. Inputting (26) into the relevant equations, the lepton masses and mixing are 
calculated as follows 

m e = 0.511 MeV, m M = 105.7 MeV, m T = 1777 MeV, 

m ni = 0.0032 eV, m„ 2 = 0.0093 eV, m m = 0.049 eV, 

s{ 2 = 0.557, 4s = 0.654, s{ 3 = 0.153, 5 l = tt, 

p 1 = -0.063 tt, (3 2 = -0.063 tt . (27) 

For a convenient comparison with the experimental data, the above results are 
converted into the commonly interested quantities in neutrino physics. They are 
such as 

Am 2 21 = 7.57 x 10~ 5 eV 2 , Am 2 2 = 2.34 x 10 ~ 3 eV 2 , 

sin 2 2B[ 2 = 0.86 , sin 2 29^ = 0.98 , sm 2 2^ 3 = 0.092 , 

m m = 0.006 eV, J l CP = , (28) 

where Am^ = m 2 a — m 2 ^, and mpp is the effective Majorana mass for neutrinoless 
double beta decay. These results in (27) and (28) are excellently in agreement with 
the recent experimental data [HE2], in particular, the value of sin 2 26[ z is identical 
with the new results of the neutrino oscillation at Daya Bay [33]. In addition, the 
valid input values in (26) are only seven, but the output values in (27) are twelve, 
the model really shows a stronger prediction power. Firstly, the lepton CP-violating 
phase angle is tt, as a result, the Jarlskog invariant J^ P is zero, i.e. the lepton CP 
violation is vanishing. Of course the reason for this arises from # 13 = in Uq. It 
should be emphasized that if a e has a imaginary part, there can be a non- vanishing 
CP violation in the lepton sector. Secondly, the two Majorana phases are the same 
but smaller. Thirdly, mpp = 0.006 eV is not too small, it is therefore promising to 
detected in the near future. Anyway, all the predicted results are expected to be 
tested in further neutrino experiments. 

In the scheme of (26), a e completely determines sin8u and sin8i3, while b e has 
only a weak effect on sm#23- The figure (3) draws the variations of the three lepton 
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Figure. 3. The graphs of the variations of the lepton mixing angles with the pa- 
rameter a e when the other parameters are fixed in (26) . 

mixing angles with the parameter a e when the other parameters are fixed in (26). 
It can be seen from the graphs that sin023 is almost unchanged, sinQyi scales down 
very slowly as a e increasing, in the meantime, sind\^ has only a weak rise toward 
right-hand side. In brief, the scenario can excellently explain that the lepton mixing 
angles are sm0 2 3 ~ 41°, sinQyi w 34° and sinQy^ w 9°. All these results essentially 
stem from the theoretical structures of the model. 

I next calculate the relic abundance of the cold dark mater neutrino v&l- The 
calculation of its annihilate cross section needs refer the two parameters tanQ g and 
y VA . By virtue of the relevant equations in (11) and (14), I respectively use the two 
mass parameters Mz 1 and M Ui to take the place of them. The Mz> value is bounded 
by both the relation Mz< ^ Zgiv^ and the experimental limits, while the value of 
M VA is determined by fitting the current observations of the relic abundance of the 
cold dark matter, namely Qh 2 = 0.112 [31]. A set of reasonable values of the two 
mass parameters are chosen as 

Mz^ = 2 TeV, M U4 = 788 GeV. (29) 

According to (11), (18) and (19), the gauge mixing angle sinQ g and the z/ 4 ^ relic 
abundance are calculated as follows 

sin0 ff = O.522, Qh 2 = 0.112. (30) 

The above results are very desired. In the figure (4), I draw the graphs of Qh 2 versus 
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Figure. 4. The graphs of Qh 2 versus M Ui for Mz> = [2, 2.5, 3] TeV, the red curve 
corresponds to the case in (29) . 



M U4 for the three values of M z > = [2, 2.5, 3] TeV, the red curve exactly corresponds 
to the case in (29). A smaller value of Mz> can not satisfy the inequality condition 
and the experimental limits, whereas a larger Mz> can lead to a very small sin9 



'9 J 



this is also unacceptable. Based on a overall consideration, a reasonable region 
of Mz> should be about 1.5 — 3.5 TeV, therefore, the corresponding values of u iL 
should lie in 700 — 1100 GeV. Because z/ 4i has no any direct interactions with the 
SM particles, however, it will be difficult for searching it in future experiments. 

Finally, I analyze the leptogenesis. The four parameters, A 4 , Ci e , c 2e , M s , have 
to be added to fulfil the calculations. The reasonable values of A4, Ci e , C2 e are all 
around 0.1. Since the role of Ci e is relatively insignificant, I directly take c<i e = C\ e 
for simplicity. The value of M5 is determined by fitting the baryon asymmetry. In 
addition, the order of magnitude of the parameter y\ is ^ 10~ 9 . In short, these 
parameters are chosen as the following values, 

y% ^ 10~ 9 , ci e = c 2e = 0.1, A 4 = 0.1, M s = 907.3 GeV. (31) 

By use of (20)-(22), the ratio of the S T decay width to the Hubble expansion rate 
as well as the baryon asymmetry are calculated as follows 

^HiT^Ms)^ ^ L9 X 10 ^ ' VB = 05 X 1(rl °- (32) 

It can be seen from the above results that the S T decay is indeed out-of-equilibrium, 
and also t]b is precisely in agreement with the current data of the baryon asymmetry 



17 




M s (GeV) 



Figure. 5. The graphs of the baryon asymmetry subjecting to Ms for the five set 
of values of the parameters [c\ e = C2 e = c e , b e ] when the other parameters are fixed 
in the context, the green curve corresponds to the case in (26) and (31) . 



[35J. Of course, Ms ~ 1 TeV is also in accordance with the previous expectation. It 
should especially be emphasized that the baryon asymmetry can still be generated 
though the lepton CP violation vanishing. The reason for this is that the fourth 
generation charged lepton and neutrino are bound to inhabit in the universe exactly 
as the model descriptions. The baryon asymmetry subjecting to Ms is drawn in the 
figure (5), in which the parameters [c\ e = ci e = c e ,b e ] are chosen as the five set of 
values, while the other parameters are fixed in the context, the green curve exactly 
corresponds to the case in (26) and (31). It can be inferred from the figure that the 
reasonable area of Ms should be about 800 — 1000 GeV. Therefore, there is a great 
chance to find the S T particle in the future. 

To summarize all kinds of the above numerical results, the model excellently 
and accurately fits all the current experimental data of the fermion masses and 
flavor mixings, and the cold dark matter and baryon asymmetry. All of the current 
measured values are correctly reproduced, while all of the non-detected values are 
finely predicted in the experimental limits. All the results are naturally produced 
without any fine tuning. In particular, the model gives a number of interesting 
predictions. A mass spectrum of all kind of the model particles is summarized as 
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Figure. 6. The diagrams of producing some non-SM particles of the model at the 
LHC, in which denotes photon. 

follows 

M v ~ 0.01 eV < M SMpartlcles ~ (0.001 - 100) GeV 

< M V4 ~ 800 GeV < (M o, M s ) ~ 1 TeV 

< {M z > , M e4 , M U4 , M di ) ~ (2 - 3) TeV < M A ~ 10 5 TeV. (33) 

In the end, I give some methods how to test the model at the LHC. On the 
basis of the model interactions, the figure (6) draws some optimum mechanisms of 
producing the model non-SM particles by the proton-proton collisions, (a) diagram 
can produce the fourth generation quarks if the colliding energy is enough effective, 
moreover, it is also a new channel of producing the SM Higgs boson H°. (b) and 
(c) respectively illustrate the productions of the heavy gauge boson Z'^ and scalar 
boson (f>°. The charged scalar boson pair S T can be produced by the (d) process. 
The fourth generation charged lepton and the cold dark matter neutrino are both 
not able to be produced in similar ways, but they can be found in the decay products 
of Zp, 0°, S T , namely Z'^ — > e^el, Z'^ — > v^pi, 0° — > ^4^4, S T — > e\v^. The loss 
of energy in the decay processes should be regarded as a definitive signal of the 
cold dark matter neutrino z/4. In particular, the (d) process produces a pair of 
S T , in succession, S T decay into e T and 1/4, the final sate leptons will eventually 
generate a matter-antimatter asymmetry. In other word, the TeV scale collider can 
experimentally produce the matter- antimatter asymmetry. In short, the model is 
quite feasible to be tested at the LHC as long as the luminance and running time 
are enough large. Undoubtedly, the best efficient methods to test the model are of 
course by the lepton-antilepton collisions at the ILC. To save space, here I do not 
go into details. A full discussion of the test of the model is planed in another paper. 

VI. Conclusions 

In the paper, I have suggested a practical and feasible particle model to account 
for the fermion flavor puzzle, the cold dark matter, and the matter-antimatter asym- 
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metry. The model has the TeV scale U(1)b-l symmetry. It contains several new 
particles with the TeV scale masses, for example, the gauge boson Z' the scalar 
bosons and S T , and the fourth generation quarks and leptons, in which the fourth 
generation neutrino is exactly the cold dark matter. In addition, the model has also 
the very heavy scalar boson A whose mass is ~ 10 5 TeV, by the new form seesaw 
A develops only the tiny VEV. This is the essential source of the tiny masses of the 
light neutrinos. These non-SM particles play key roles in the model. The model 
can not only excellently explain the fermion masses and flavor mixings, but also ele- 
gantly accommodate the cold dark matter and leptogenesis at TeV scale, moreover, 
it gives some interesting predictions. The theory can perfectly integrate three party 
of the flavor physics, the cold dark matter and the matter-antimatter asymmetry, 
therefore, it is quite deserved to be tested in future experiments on the ground and 
in the sky. Finally, I believe that the new non-SM particles, including the cold dark 
matter neutrino z/4, are possible to be discovered at the LHC and ILC in the future. 
In a word, all kinds of the experiments toward these goals will not only provide us 
more information about particle physics, but also help us understand the mysteries 
of the universe. 
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